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and Richard P. Shaw 

SUMMARY 


Minimum critical Reynolds numbers for laminar instability of two- 
dimensional stagnation flows over a wall cooled by normal fluid injection 
are determined theoretically. The individual as well as the net effects 
of wall temperature and normal fluid injection on the stability charac- 
teristics of the laminar boundary layer are thereby investigated. In 
order to study the effect of the favorable pressure gradient in this type 
of flow, comparison is made between results of calculations obtained 
herein and results of similar calculations for the flow over a flat plate. 
A comparison of various stability criteria is included, and, for certain 
cases of normal mass-flow injection, a slight appropriate extension in 
the approximate stability criteria is made. 

A particularly interesting result obtained here is that, at fixed 
zero or low rates of fluid injection, the stabilizing effect of cooling 
of the wall is much greater for the flow over a flat plate than for the 
stagnation flows. In fact, over a highly cooled wall with zero or low 
rates of injection, the stagnation flows were found to have a minima mi 
critical Reynolds number based on momentum thickness lower than that of 
the corresponding flat-plate flows. 


INTRODUCTION 


The purpose of this study is to investigate theoretically the sta- 
bility of the laminar boundary layer in flow with an axial pressure gra- 
dient over a wall cooled by normal fluid injection. In particular, the 
separate as well as the simultaneous or net effects of wall temperature, 
of a small normal mass flow at the wall, and of a pressure gradient on 
the stability characteristics are determined. Such an Investigation is 
of practical Interest in connection with, for example, the transpiration, 
or sweat cooling, of rocket walls, turbine blades, and airfoils. A study 
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of the stability of the laminar boundary layer is of practical interest 
since, at least under conditions of a low free-stream turbulence, insta- 
bility of the laminar boundary layer appears to be a necessary condition ' 

for transition to turbulence. The main advantages of maintaining a lami- 
nar boundary layer under these conditions would be a relatively low skin 
friction (provided that separation does not occur, cf . ref. 1) and a 
relatively low required normal mass flow of the coolant to maintain a 
given wall temperature. 

The simultaneous effect of cooling of the wall and of a small nor- 
mal mass flow at the wall is of especial interest, since the separate 
effects tend to oppose each other. In particular, lowering the wall tem- 
perature tends, by itself, to stabilize the laminar boundary layer, while 
a normal mass flow at the wall has,' by itself, a destabilizing effect. 

It is therefore significant to determine under which conditions and to 
what extent the net effect of transpiration cooling is stabilizing or 
destabilizing. 

The stability of the laminar boundary layer has already been studied 
to an appreciable extent in the literature, but most of the studies have 
been carried out for an impermeable wall, that is, for zero normal mass 
flow at the wall. A survey (as of 1952) of theoretical and experimental . 

results related to such studies can be found in reference 2, while theo- 
retical results of additional and later studies are summarized in refer- 
ence 3-. The earliest stability investigations of flow with normal mass ' 

flow at the wall were concerned with the effects of suction. References 4 
to 6 show the stabilizing effect of suction in increasing the minimum 
critical Reynolds number of the incompressible laminar boundary layer over 
a flat plate with uniform suction and with a suction velocity inversely 
proportional to the square root of the distance from the leading edge. 

The laminar-boundary- layer solutions -used for this purpose were the exact 
or numerical solutions obtained In references 7 and 8. 

Investigations of the stability of flows with normal fluid injection 
have been made in references 9 to- 11. In reference 9 , the stability of 
the laminar boundary layer over a flat plate with a transpiration-cooled 
wall maintained at a uniform temperature is studied. Although the flow 
is here essentially at low speeds (zero or negligible Mach number), den- 
sity and temperature variations in the flow are taken into account. The 
Injection velocity at the wall is assumed to vary inversely as the square 
root of the distance from the leading edge, and it is shown that by a sim- 
ple transformation the incompressible solutions of reference 7 can be 
applied. Results for the minimum critical Reynolds number show explicitly 
the stabilizing effect of cooling of the wall and the destabilizing effect 
of normal fluid injection. In reference 10, the compressible flow over a 
flat plate, with nonzero Mach number and uniform mass flow (both suction 
and injection) at the wall, is considered. The boundary-layer solutions 
used are, however, obtained by means of the Karman-Pohlhausen method in *v 
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conjunction with, fourth-degree velocity profiles and may, therefore, 
yield appreciable quantitative errors in the stability calculations (cf., 
e.g., ref. 12). Reference 11, which is quite recent, presents and employs 
exact flat-plate boundary-layer solutions which are essentially an exten- 
sion to high-speed flows, that is, to nonzero Mach numbers, of the solu- 
tions used in reference 9* The use of liquid air and water as well as of 
air as coolants is studied. The stability calculations in reference 11 
consist in determining the (uniform) wall temperature as a function of 
Mach number and of the injection mass-flow parameter required to stabilize 
completely the laminar boundary layer. These results may be considered as 
extensions to flow over a flat plate with normal fluid injection of the 
corresponding results of references 13 and 14 for an impermeable flat 
plate. The destabilizing effect of the normal fluid injection is ill dr. 
trated here by its effect in lowering, for a given (supersonic) Mach num- 
ber, the wall temperature required to stabilize the flow completely. Ref- 
erence 15 contains results of stability calculations for incompressible 
flow in a pressure gradient with suction or injection, based on Schlichting* 
one-parameter integral method . The emphasis here is on the effect of a 
change in boundary-layer thickness on the stability of the lamina r boundary 
layer. 


In all of the foregoing research, except reference 15, on stability 
of laminar flows with normal mass flow at the wall, only the case of a 
zero pressure gradient (such as flow over a flat plate) has been consid- 
ered. In the present Investigation, the effect of a pressure gradient is 
taken into account by treating two-dimensional flows In the vicinity of a 
forward stagnation point over a transpiration-cooled wall. Such flows are 
not only of theoretical interest, but are also of practical interest, 
since they occur, for example, near the leading edge of a blunt-nosed 
object in subsonic flow. Moreover, accurate boundary-layer solutions for 
such flows have been derived and are tabulated in references 1 6 and 17 . 
These tabulations are particularly suitable for the stability calculations 
to be performed here, since they include values of the first and second 
derivatives of the velocity and temperature profiles. The flows consid- 
ered here are low-speed or zero Mach number flews, but they are neverthe- 
less compressible in the sense that density and temperature variations are 
taken into account. In order to determine properly the effect of the 
favorable pressure gradient inherent In the stagnation flows, the stability 
calculations have also been performed for the flow over a flat plate, with 
essentially the same values of wall-temperature and Injection-mass -flow 
parameters. Tie corresponding laminar-boundary- layer solutions for flow 
over a flat plate are also tabulated in reference 16. By calc ula ting the 
minimum critical Reynolds numbers as functions of the wall-temperature 
and injection-mass-flow parameters for both stagnation and zero-pressure- 
gradient flows, the separate and simultaneous effects of wall temperature, 
injection mass flow, and pressure gradient are studied in the present 
investigation. 
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The stability criteria to be used in this study are primarily the 
approximate two-dimensional criteria developed by Lin and Lees (refs. 18 
to 20) with modifications along the lines of Bloom (refs. 14, 21, and 
especially 22) . As will be seen, however, a slight extension in these 
criteria is necessary for certain types of cases, including the flow 
over a flat plate with a normal injection mass flow at the wall. This 
extension is made herein. Subsequent to the completion of most of the 
calculations, a report by Dunn and Lin (ref. 23) appeared, in -which a 
slightly modified set of stability criteria was developed. This report, 
moreover, included a treatment of three-dimensional disturbances. In the 
present investigation, consequently, slightly more than half of the cases 
were recalculated accordingly, and the results were compared with those 
based on the previous criteria. The comparison indicates that, although 
quantitative differences appear, the conclusions for the flows investi- 
gated herein based on either set of criteria remain quite similar. It 
will be seen, moreover, that for the low-speed though compressible flows 
treated herein, it is sufficient to determine the stability character- 
istics for only two-dimensional disturbances, since the flows will be 
more stable with respect to three-dimensional disturbances. 

The present report Is divided into three sections. In the first 
section the stability criteria to be used in this investigation are dis- 
cussed explicitly. The mean-flow boundary- layer solutions Investigated 
herein and the application of the stability criteria to them are then 
discussed in detail in the second section. Finally, the third section 
presents the results of the stability calculations and the conclusions 
drawn from them. 

This research was carried out at the Polytechnic Institute of 
Brooklyn Aeronautical Laboratories under the sponsorship and with the 
financial assistance of the National Advisory Committee for Aeronautics. 
The authors hereby express their thanks to Professor M. Bloom for his 
helpful discussions. 


SYMBOLS 


c 


C P 


dimensionless complex phase velocity of disturbance, 
c r + icj_ (c r and Cj_ real) j for a neutral disturbance, 
c = c r 

value of c at minimum critical Reynolds number 
specific heat at constant pressure 


f 


function of tj related to velocity profiles by equation (Bl) 
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f w value of f at ■wall, proportional to normal mass flow at 

wall (negative for injection) (cf. eqs. ( 29 ) to ( 32 )) 

G(z,A) function in stability analysis, defined by equation (j) 

G min (A) minimum value of G(z,A) with, respect to z for a fixed 
value of A 

J(c) function in stability analysis, given by equation (5) 

j(z,A) = l/G(z,A) 

J miT1 (c) = l/G^-t n (A) (with A as a function of c in accordance with 

eqs. (8) or (10)) 

K factor in expression for momentum thickness (eq. (35)) 

k coefficient of heat conductivity 

L boundary- layer length (such as momentum, thickness) used in 

stability calculations 

Z characteristic main-stream length (such as chord) 

M Mach number 

Mi = M 1 (a/cx 1 ) 

m exponent in equation ( 25 ) 

Pr Frandtl number 

R critical Reynolds number, u^L/v^j that is, values of uqL/vq 

along neutral stability curve 

R 1 ° R ( ct '/ ct l) 


R b 




Reynolds number based on conditions at reference point out- 
side of boundary layer, u^l/v^ 

Reynolds number based on momentum thickness, u qE>i/ v i 


T absolute temperature 

coolant temperature 
t time 


u 


velocity in x-direction 
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v w normal velocity component at -wall 

v = u/u-l 


x,y 

z 

z 


a 


coordinates parallel and normal to surface, respectively 
coordinate perpendicular to x and y 

variable used in stability analysis (see eqs. (4) and (ll)) 

value of z corresponding to minimum value of G(z,A) for 
a fixed value of 

dimensionless wave number of two-dimensional disturbance 

(eq. (1)) 



p proportional to dimensionless wave number in three-dimensional 

disturbance (eq. (22)) 

T)T V angles defined in figure 6 


e 


n 

$ 

a, 

K 


momentum thickness of boundary layer 

e 

exponent in heat-conductivity— temperature relation, k « T 
variable defined by equation (26) 
temperature-difference ratio, (I - vJlfa - T„) 

K 

exponent in specific-heat — temperature relation, Cp « T 


stability profile parameter (see eqs. (8) and (10)) 
jjt coefficient of viscosity 

v kinematic viscosity, p/p 

| dimensionless distance along wall, x/Z 

mass density 


P 
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4 r (z)^i(z) 

9(1) 


real and imaginary parts , respectively, of function 
related to Tietjens function 

injection mass-flow parameter, ^p^v^/p^u^ \Jb^ 


qj average value of cp(g) over length l 

to exponent in viscosity-temperature relation, p « 

Subscripts: 


1 values at local outer edge of boundary layer 

b values at reference point outside boundary layer, namely, 

at x = 2 

c values at u/u^ = c 

c Q values at u/up = c Q 

min minimum critical value, for example, R_p n 

w values at wall 

A prime in the main text denotes differentiation with respect to 
y/L (except for ; in appendix A, a prime denotes differentiation 

with respect to yj in appendix B, it denotes differentiation with respect 
to tj . 

A bar denotes values based on equations (10) to (12) for modified 
two-dimensional stability criteria (except for qj) . 


STABILITY - CRITERIA 


To aid in the physical interpretation of the results to be obtained 
herein, the stability criteria upon which these results axe based will 
first be briefly reviewed and discussed. Such a brief review appears 
particularly worthwhile in view of the recent literature on this subject. 


Curve of Neutral Disturbances 

The stability criteria to be applied in the present investigation 
are based primarily on those developed in references 19 and 20j these 
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stability criteria are obtained by an extension to compressible flew of 
the two-dimensional theory developed in reference 18 for incompressible 
flow. The stability criteria are based on the amplification or decay 
with time of small disturbances superimposed on the steady-state boundary 
layer, which is assumed to be known. The disturbances Q considered are 
of the form 


Q(x,y,t) = F(y) e i ( a A')( x ~ cu l t ) (1) 

(The length L and the velocity u^ are introduced in eq. (l) so that 

a and c as defined there are nondimens ional. } The wave number a is 
real and positive and is inversely proportional to the wavelength of the 
disturbance. The disturbance velocity c may be complex. Thus, 
c = c r + icj_, where c r Uj_ is the phase velocity of the disturbance, 

while c-j_ determines the time rate of amplification (c^ > 0) or decay 
(ci < 0) of the disturbances. The case Cj_ = 0 defines a neutral dis- 
turbance and determines the bounds between stable and unstable flows. 

By applying the Navier-Stokes equations, in conjunction with approxi- 
mations pertaining to boundary- layer flows, it is found that the stability 
of a particular boundary layer at any given point x along the flow 
depends only on the local properties (such as the velocity and temperature 
profiles) of the boundary layer at this point and not, for example, on the 
previous history of the flow. This has been found to be true, as a first 
approximation, even when axial gradients of the disturbances and the nor- 
mal velocity component of the steady-state flow are taken into account in 

the disturbance equations (refs. 2k to 26). ^ Thus a set of linear homo- 
geneous partial differential equations in the disturbances is obtained, 
with solutions in the form of equation (l) . When equation (l) is inserted 
into the partial differential equations, a set of linear homogeneous 
ordinary differential equations, with y as the independent variable, is 
obtained. The condition for the existence Of a nontrivial solution, under 
appropriate homogeneous boundary conditions, then yields a relation 
between a and the eigenvalues c, with the Reynolds number R as a 
parameter. This relation can be written as two equations of the form 


c r = c r (a,R) 

Cj_ = Ci(a,R)^ 


( 2 ) 


^It may be noted in this connection, however, that according to 
reference 26 the normal-velocity terms in the basic differential equa- 
tions are of a larger order of magnitude than the axial-gradient terms. 
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By setting = 0, which, is the condition of neutral stability, a rela- 
tion between a and R, of which a typical form is shown in figure 1, 
is obtained. For a supersonic main stream (Mi > 1) at the edge of the 
bou nd a r y layer it is assumed that the disturbances mus t be subsonic, 
which requires, in such cases, that 

c r > 1 - (1/Mi) (3) 

The lower branch of the neutral curve of a versus R is asymptotic 
to the R-axis , that is, R-*co as cc -» 0. The asymptotic behavior of 
the upper branch depends in general on the vanishing or nonvanishing of 
the quantity (S/^y)(p Bu/By) in the local boundary layer. Thus, the 
upper branch also approaches the R-axis as R -» « provided that 
(B/By) (p Bu/By) / 0 in the boundary layer j otherwise, if 
(S/By) (p Bu/By) = 0 at u/ui = c s say, then the upper branch has a 
positive asymptote, namely a -»<%, where <xg corresponds to c = c s 
in the characteristic equation for an infinite value of R. The portion 
of the curve inside the two branches has been shown to correspond to 
pairs of values of a and R for which the boundary layer will be 
unstable, while for the pairs of values of a and R outside the loop, 
* the flow will be stable. 


Minimum Critical Reynolds Number 

Of particular interest in the present investigation is the minimum 
value of R, called the minimum critical Reynolds number and denoted as 
Rmin j on the neutral curve. From figure 1 it can be seen that R^n is 
the Reynolds number below which all (small and subsonic) disturbances, 
regardless of wavelength, will be damped out. Thus, a necessary condi- 
tion for unstable boundary-layer flow, at least with respect to fiman 
two-dimensional disturbances of certain wavelengths, is that the actual 
Reynolds number at same station exceed R^ji at that station. If the 
free-stream turbulence is sufficiently low, then this appears to be a 
necessary condition for eventual transition from a laminar to a turbulent 
boundary layer. Actual transition, however, occurs downstream of the 
point where the minimum, critical Reynolds number is first exceeded and 
depends on additional factors such as the actual size of the disturbances 
and hence on their rates of amplification as well as on nonlinear effects. 
Even in this respect, however, the value of n may have, special sig- 
nificance, for it is stated in references 9 and 20 that the maximum rate 
of amplification of the selfexcited boundary- layer disturbances varies 

roughly as l/yj^dj) min' Thus, the magnitude of the minimum critical 

Reynolds number can serve at least as a qualitative measure of the degree 
of stability, and hence of the tendency for transition, of a given boundary- 
layer flow. Since the present analysis is concerned pr imar ily with the 
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effect of various parameters, such, as the wall temperature and normal 
injection mass flow, on the stability of the laminar boundary layer, a 
convenient’ and suitable quantitative means of carrying out such an anal- 
ysis is to calculate and compare the minimum critical Reynolds numbers 
for the various cases. 


Two-Dimensional Stability Criteria and Calculation 
of Neutral Curve 

According to the development in references 19 and 20, the neutral 
curve a = a(R) involves the variable z defined as 


» - (y/L) 

V V =/ V 1 / 


(M 


where the subscript c denotes values at the point where u|u-^ = c, 

called the inner critical layer. The prime here denotes differentiation 
with respect to y/L. The actual location of the critical layer depends 
largely on the function (using the notation of ref. 22) 


J(c) 


“it 


(u/ui) w ' C (T c /I 1 ) a 

[(“/"■!) '] 5 (V"1> 



1 


(u/Ul) v '=(T <: /T 1 ) 

'(u/ai)- 

(ipiY 

[(V"i)'] c Wi> 

S u / u 1) ' 

(T/Ti) _ 


^observing that p/p^ = T p/f )• Along the neutral curve, 

J(c) = j(z,A) 


(5) 


(6) 


where j(z,A) is a certain function of z and a boundary- layer profile 
par ame ter A defined below by equation (8) (functions J(z,A) and J(c) 
correspond, respectively, to' Tjf ± (z,c) and v (c) in ref. 23). The func- 
tion J(z,A) can be expressed as (cf. refs. 22 and 23) 


1 

j(z,A) 


=G(z,A) 




(*r 2 4- tf) 

a + M 


(?) 
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where $ r (z) and ®j_(z) are the real and ima ginary parts of a function 

related to the Tietjens function and can be found tabulated or plotted in 
references 18, 19, 20, or 23* The parameter A is given by 



Substituting for (y/L) c according to equation ( 8 ) into equation (4) and 
solving for aR yields 


( u / u i ) c ' c5 ( 1 + > 0 5 


An accurate determination of the wave number a. according, for example, 
to reference 20 (eq. ( 27 ) and appendix) involves the calculation of inte- 
grals with certain types of singularities (cf . also ref. 23 , where this 
procedure is somewhat simplified) . 

For any given boundary-layer flow, equations ( 6 ) and ( 9 ) and an equa- 
tion for a determine points on the neutral curve in the aR plane. 

Thus, for a chosen value of c, j(c) and A(c) follow from equations ( 5 ) 
and ( 8 ), respectively, while j(z) follows from equation ( 7 ). Values of 
z satisfying equation ( 6 ) are then determined. If there are no values 
of z satisfying equation ( 6 ), then no neutral disturbance with the par- 
ticular chosen (dimensionless) phase velocity c will exist for the given 
flow. If values of z satisfying equation ( 6 ) for a given value of c 
do exist, then it will be found in general that two such values occur. 

This is due primarily to the nature of the G(z,A) curves (cf. fig. 2). 
For each of these values of z, aR can be calculated from equation (9) , 
while a can be determined from an additional equation. This yields a 
pair of values of a and R. 


In a recent investigation, including three-d ime nsional disturbances, 
Dunn and T in (ref. 23) have, for greater accuracy, slightly modified the 
above two-dimensional stability criteria. The parameter A (to be denoted 
here by a bar) is now expressed as 



Moreover, from a modified definition of z, aR (also denoted by a bar) 
is now given by 
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With the use of equation (10), equation (ll) can also be -written as 

_ . (12) 

c3(l + A) 2 

Although equations (10) and (12) may in certain cases yield substantially 
different quantitative results from those obtained by equations (8) and 
(9) , the new equations are on the whole essentially similar to equa- 
tions (8) and (9 ) > and in cases of very moderate temperature changes in 
the boundary layer the two sets of stability criteria will even yield 
similar quantitative results . A more detailed comparison of these two 
sets of criteria is given in subsequent sections. 

The calculation of the neutral curve as described herein can be 
lengthy and tediqus. If, however, it is desired only to calculate the 
minimum critical Reynolds number, corresponding therefore to only one 
point on the neutral curve (cf. fig. 1), then, with certain approxima- 
tions, the calculations can be considerably simplified. Such approxima- 
tions are particularly useful if, as in the present investigation, it is 
desired to determine the effect of various parameters, Bince the consist- 
ent use of an approximate set of stability criteria should suffice to 
give a quantitative as well as a qualitative indication of these effects. 

The minimum critical Reynolds number may be assumed to correspond 
to that value of the wave velocity c for which there will be a single, 
and only a single, value of z satisfying equation (6) (refs. 22 and 25). 
This is equivalent to assuming that R m j_ n corresponds to that value of z 

for which G(z,?\) will be a minimum for a fixed value of ~h (cf. fig. 2). 
This assumption appears reasonable since it has been seen (cf . fig. 2) that 
the minimum point on a G(z,7\) curve marks the boundary below which neu- 
tral disturbances cannot exist. It may also be noted that the value of c 
corresponding to a minimum value of G will usually be found to be the 
maximum value of c for which equation (6) can be satisfied. This appears 
to be due primarily to the fact that J(c) (s> l/G(c)), at least in the 
vicinity of its intersection with the G m j_ n (c) or J^c) curve, 
increases monotonically with c. These arguments are based on neglecting 
the variation of A(c) with c. This, however, appears to be justified, 
since, at least in the vicinity of the critical layer, it has been found 
that 7\(c) varies relatively very slowly with c. 
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In cases in which |A[ « 1, it follows from equation (7) that 
G(z,A) will be a minimum for a given value of c (or A) at approxi- 
mately the value of z, namely z = 3.22, at which $^(z) is a maximum . 

Since at ' z = 3-22, = 0.580 and 4> r = 1.50, it follows from equa- 

tion ( 7 ) that for small values of | A | 

<WM ~ cr (13 > 

Equation (6) then becomes 


(1 - 2A)J(c) = 0.580 (l4) 

Equation (l4) is the familiar approximate equation which has been devel- 
oped and used in references 18 and 20 and has been applied in most approxi- 
mate stability calc ulat ions. 

In certain of the cases studied in the present investigation, namely, 
the flow over a flat plate with appreciable normal fluid injection at the 
wall, large negative values of A (of the order of magnitude of -0.7) 
are obtained, and equation (13) (and, hence, also eq. (l4) ) is then no 
longer strictly valid. In place of equation (13) , minimum values of 
G(z,A) with respect to z for various fixed values of A can be calcu- 
lated straightforwardly by means of equation (7). This has already been 
carried out by Bloom (ref. 22) for large positive^ values of A. For 
negative values of A, G(z,A) is shown plotted in figure 3 as a function 
of z for various fixed values of A when A > -1. From figure 3j for 

each fixed value of A, the values of z (to be denoted as z D ) and of 

G (denoted as Gm-tn ) at which G is a minimum are determined. The values 

of z Q and G m -; n as functions of A are given in figures 4 and 5> 

respectively. It may be noted from figure 4 that z Q varies fairly 
slowly (from z Q = 3*10 at A = 0.1 to z D = 3*74 at A = -0.8) with 
A in the range shown; G Tn -r ri (A), on the other hand, varies rapidly at 

appreciable negative values of A and becomes infinitely large as A -> -1 
(fig. 5)- For very low values of J A | , it is seen from figure 5 that 
(^^(A) is given approximately by equation ( 13 ) • 


In carrying out calculations for flow over an impermeable flat plate 
at high Mach numbers and large rate of cooling on the basis of equation (8), 
very large positive values of A were obtained by Bloom. Using equa- 
tion (l0), however, Dunn and Lin (ref. 23) subsequently obtained small 
(positive) values of A for all such cases of flow over an impermeable 
flat plate . 
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Physical or Geometric Significance of A 


It may he of interest to note that the stability parameter 'A has 
a fairly simple physical or geometric significance. First, however, it 
should he observed that equation (10), which is the more accurate equa- 
tion for A, is, at least for cases of moderate temperature changes within 
the boundary layer, qualitatively quite similar to equation (8) . It can, 
in fact, he shown (appendix A) that if v is assumed constant, then an 
expansion of the velocity profiles about the critical layer u/u-j_ = c 

indicates that to a first approximation both A and X will be propor- 
tional to the second derivatives of the velocity profiles and will have 
the same order of .magnitude. It will subsequently be seen that _f or all 
of the (low-speed) caseB treated herein the values of A and A are 
essentially similar. (For cases of high Mach .number and very large tem- 
perature changes within the boundary layer, however, eq. (10) may yield 
significantly different results from eq. (8) because of the presence of 
the kinematic viscosity v terms in eq. (10)). 


It is seen from the expansion about the critical layer (appendix A) 
that A (or *h) depends on the curvature of the velocity profiles. The 
geometric significance of A can be seen particularly clearly from the 
approximate equation (8). From figure 6, in fact, it can be readily 
seen that, according to equation (8), if the velocity profile Is plotted 
(in the customary manner) with y as ordinate and u as abscissa, then 
at any point on the profile 


A = jaa - 7 -- - l 
tan 7 


'w 


(15) 


where tan y^ is the slope of the velocity profile (dy/du) = l/(du/dy) 

at the wall, while tan y is the slope of the secant from the origin to 
the given point on the velocity profile. Thus, 1 A I will be small if the 
velocity profile is almost linear near the wall; however, A will, for 
example, have a large negative value if y is appreciably smaller than 
7 W and hence if the velocity profile has ah appreciable convex curvature 
(cf. fig. 6). The velocity profiles for flow over a flat plate with 
fluid infection at the wall are found to be of the latter type. On the 
other hand, it will be seen that the velocity profiles for stagnation 
flows with or without fluid injection all have only slight curvature near 
the wall and hence are associated with small values of | A I . 


Calculation of c 0 

The minimum critical Reynolds number depends strongly on the loca- 
tion of the inner critical layer, as given by the value of c. In order 
to determine c Q , the value of c corresponding to the minimum critical 



MCA TN 4037 


15 


Reynolds number , J(c) 
c (or equivalently of 
or (10) , respectively, 
values of A, (e) 


and A(e) are first calculated as functions of 
(y/L) c ) in. accordance with equations (5) and (8) 
For each value of c with the corresponding 
is obtained from figure 5- Ihe value of c Q is 


the value of c for which 


J(c) 


1 


( 16 ) 


Equation (l6) can be solved by plotting J(c) and J TO in (c) and obtaining 
their point of intersection. Examples of such curves are given in fig- 
ure 7.3 This method of determining c Q , in conjunction with figure 3 > , 

has been consistently applied in the present report for all of the cases 
studied herein, (in refs. 9 to 11, where the stability of the laminar 
boundary layer over a flat plate with normal fluid injection is treated, 
no mention is made of the actual numerical values of A obtained therein.) 


Calculation of Minimum Critical Reynolds Number 

In reference 9, the following approximation for the wave number a 
has been used for the calculation of the minimum critical Reynolds number 
of the laminar boundary layer over a sweat -cooled flat plate: 


a 


( u / u l)w' 
1.2 



(17) 


This expression has evidently been obtained from equation ( 27 ) of refer- 
ence 20 by neglecting all terms proportional to a (or powers of a) 
on the right side and replacing the quantity u - L by the constant 1.2. 
The latter approximation for u - L may be adequate at least for low 
speeds (Mj_ « 0 ). The use of a simple approximate expression such as 

equation ( 17 ) for a introduces a considerable simplification in the sta- 
bility calculations. In cases for which the wall temperature ratio T w /Tj_ 
does not differ greatly frcam unity, equation ( 17 ) is substantially si m ila r 
to the more frequently used approximation developed in reference 20, 
namely. 


a = ( u / u l) * c o Z 1 " “ c o) 2 ( l8 ) 


3 

In the actual calculations, these curves were all plotted on more 
suitable individual scales . The chief purpose of figure 7 is for the 
subsequent comparative discussion of the results . 
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For M]_ = 0, equation (l8) reduces to the approximation for a developed 
in reference 18 for incompressible flow. 


By inserting expression (17) for a into equation (9), the following 
expression is obtained for the minimum critical Reynolds number: 


1 . 2 Zo :i (v e c,/v 1 )[(u/ui) ¥ ’] 2 (I w /Ii)(l - Co ) 2 

°oV/"l) c 0 ’( 1 + *c 0 ) 3 [ 1 - - =o)f /2 


Similarly, equation (12) in conjunction with equation (17) yields 


c 0 f 

R mln — r y/2 ^ 20) 

+ S =o) t 1 - *1 4 - °°) 2 ] 7 

In equations (19) and (20) the subscript c Q denotes values at u «= c Q . 

Equation (20) is, in general, essentially similar to equation (19) , espe- 
cially in cases of small thermal changes in the boundary layer (for which 
v c j v 2 _ ~ v w/ v l) and of small curvature of the velocity profiles (for which 

| A} or (a| « 1, while (u/u-^) * ** ( u / u l) w ')*^ ^ should be noted, how- 

ever, that the value of c Q to be used in equation (20) will not neces- 
sarily be exactly equal to that in equation (19) , since the former should 
be based on expression (10) for X, while the latter should be based on 
equation (8) . 

It may be noted that in cases of small values of | A ( and c Q , with 
T w /Tj_ close to unity, equation (19) or equation (20) reduces (except for 

a constant factor of 1.2 X (5*2)^ = 59*5; instead of 25) to the following 
relatively familiar approximate expression developed In references 18 
and 20: 


*Wn 


g ( v ° 0 / v i)( u K) v ' 

=^[1 - M^CL - c 0 ) 2 ] 


( 21 ) 


^Even In cases of fairly large curvature, equations (19) and (20) 
tend to be similar even quantitatively. For example, if A(« A) is nega- 
tive, then the additional (l + A Cq ) factor in the denominator of equa- 
tion (19) tends to make Rmi n greater than R^n, but the fact that in 
such a case ( u / u i) ’ > ( u / u i) ' tends according to equations (19) 

and (20) to compensate (sometimes even overccmpensate) for this effect. 
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In the present invest igation, all of the calculations were first carried 
out on the basis of equation ( 19 ) in conjunction with equation (8), and 
subsequently, for comparison, a majority of the cases were recalculated 
on the basis of equation (20) in conjunction with equation (10) . In all 
of these cases, the -value of was M]_ = 0. 


Three-Dimensional Disturbances - Extension of Squire’s Theorem 

The stability criteria discussed in the foregoing sections are based 
on two-dimensional disturbances, represented mathematically by equations 
of the form of equation (l) . Actual disturbances may, however, be three 
dimensional. Such disturbances, for boundary- layer flows, can be repre- 
sented in the form 


d(x, y ,z,t) ,y (y)e i[WL)»(eA)z-(»AK=t] (22) 

where p (as well as a) is real and positive. 

For incompressible. flow, it has been shown by Squire (ref. 27) that 
the stability equations for three-dimensional disturbances in the boundary 
layer can, by single suitable transformations, be reduced to the equations 
of two-dimensional disturbances with wave number a,^ and Reynolds num- 
ber Rp defined by 


2 2 , □ 2 

On = a, + p 

Rj_ = R(a/aq) 


(23) 


Since Rp < R, it is concluded that any instability which may be present 

for a two-dimensional disturbance will also be present for a three- 
dimensional disturbance at a higher Reynolds number . Consequently, incom- 
pressible flows are more stable with respect to three-dimensional disturb- 
ances than with respect to two-dimensional disturbances, and it is 
therefore sufficient, in practice, to analyze the stability of incompres- 
sible flows with respect to only two-dimensional disturbances. 

For compressible flow (at Mach numbers below the hypersonic range), 
it has been recently shown (refs. 23 , 28 , and 29 ) that by transformations 
essentially analogous to those for incompressible flow, the stability 
equations for three-dimensional disturbances in the boundary layer can, 
once again, be reduced to those for two-dimensional disturbances with wave 
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number ctq_ and Reynolds number Rq as given by equations (25) How- 
ever, the Mach number Mq corresponding to these quasi-two -dimensional 
disturbances will be 


&1 = Mq(a/ct,q) (24) 

Because of the transformation of the local Mach number as given by equa- 
tion (24) , the quasi-two-dimensional disturbances do not represent proper 
two-dimensional disturbances, and Squire's theorem on the greater sta- 
bility of a given flow with respect to three-dimensional disturbances 
than with respect to two-dimensional disturbances is no longer strictly 
valid (cf. ref. 25). On the other hand, for compressible flows at zero 
Mach number (i.e., at low speeds), such as those to be investigated in 
the present report, there is no transformation of Mach number, and 

£ 

Squire's theorem remains valid. Thus, for the cases studied herein, 
it is sufficient to analyze the stability with respect to only two- 
dimensional disturbances. 


MEAN-FLOW OR STEADY-STATE SOLUTIONS 


The stability criteria described in the section entitled "Stability 
Criteria" depend on the particular steady-state or mean-flow boundary- 
layer solution considered. It has been Been, moreover, that the stability 
criteria are fairly sensitive to first and especially to second deriva- 
tives of the velocity profiles (cf. eqs. (5))* Consequently, it is neces- 
sary that the mean-flow solutions considered be sufficiently accurate with 
respect to these derivatives. In the present study, such solutions are 
required for the case of compressible (i.e., variable density and tempera- 
ture) flow with heat transfer In a pressure gradient with a small normal 

- ■ “ 

The energy equation, not considered in reference 27 for incompres'- 
sible flow, appeared at first (ref. 29 ) to yield additional terms which 
do not permit a strict reduction to quasi-two-dimensional-disturbance 
equations. However, it is shown in reference 25 that these terms will, 
at least for moderate Mach numbers, be as small as other terms which have 
been neglected. 

^This can be further verified from the condition (9«12) of refer- 
ence 25 under which compressible three-dimensional disturbances have 
a minimum critical Reynolds number ' higher than that of two-dimensional 
ones. If Mq = 0, then this condition will be seen to be identically 
satisfied. Moreover, at nonzero but subsonic free-stream Mach numbers, 
three-dimensional disturbances, according to reference 25, will usually 
have a minimum critical Reynolds number higher than that of two- 
dimensional ones, except possibly under conditions of extreme surface 
cooling. 
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injection velocity at the wall. Only a limited number of sufficiently 
accurate solutions for such fairly general cases appear to have been 
obtained thus far. A class of boundary-layer solutions of particular 
interest is that for flow over a transpiration-cooled wall with a uniform 
coolant, as well as a uniform wall, temperature. Approximate solutions 
for such flows, with arbitrarily prescribed pressure gradients and Mach 
numbers (below the hypersonic range), have been developed in reference 3 0 . 
These, however, are based on the assumption of fourth-degree profiles and 
may therefore not be sufficiently accurate for stability calculations 
(cf. ref. 12), except, perhaps, for the stagnation flows explicitly cal- 
culated therein (cf. ref. 31). Solutions of greater accuracy for an 
arbitrarily prescribed pressure gradient can be developed without dif- 
ficulty from the analysis in reference 1. It would merely be necessary, 
for this purpose, to apply sixth-degree velocity profiles (emitting the 
boundary condition (l7d) in ref. l) instead of the seventh-degree velocity 
profiles used therein. Another set of solutions for flow over a sweat- 
cooled wall which have been developed is that based on low-speed, but 
compressible, flow with the local velocity outside the boundary layer 
proportional to a power of the downstream distance along the wall (refs. 16 , 
17, and 52). Since, as will be seen, these solutions are particularly 
suitable for the purposes of the present investigation, these are the 
solutions the stability of which will be analyzed here. 


Boundary-Layer Solutions for "Wedge Flows" 

Flows for which the local velocity at the outer edge of the boundary 
layer is proportional to a power of the downstream distance along the wall 
have been called wedge flows since, in incompressible flow, they can repre- 
sent the potential flow over wedges of various angles (see, e.g., ref. 35) • 
Such flows can be characterized by the following nondimens ional equation 
for the velocity distribution outside the boundary layer: 


“a7“b - i m (=5) 

where m is a given constant, and the subscript b denotes (reference) 
conditions at | = 1 or x = 2 outside the boundary layer. The pres- 
sure gradient in such a flow will be positive (adverse) or negative 
(favorable) according to whether m < 0 or m > 0. Such flows were 
originally investigated for incompressible flow over an impermeable wall 
in reference 34. The case m = 0 represents the flow over a flat plate 
(zero pressure gradient), while the case m = 1 represents the subsonic 
flow in the immediate vicinity of a two-dimensional forward stagnation 
point. These cases, therefore, have particular physical significance. 
For zero Mach number, as will be assumed here. 


T l/^b = Pl/Pb = 1 


along the entire flow 
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In reference 32 , boundary-layer solutions for flows characterized 
by equation (25) (especial l y with m = 1) over a sweat-cooled wall were 
developed for incompressible flow with constant fluid properties. Such 
solutions were subsequently extended in references 35 > 16, and 17 to low- 
speed compressible flows with variable fluid properties. The wall tem- 
perature in all of these flows is assumed to be uniform. The mathematical 
advantage of treating such flows is that, by a fairly simple transforma- 
tion, the partial differential equations of the boundary layer can be con- 
verted exactly into ordinary differential equations. Thus, a class of 
similar solutions is obtained, with the velocity u/uq and tempera- 
ture T/Tq profiles as functions of a single variable q defined by 


4 



(26) 


In references 16 and 17, the boundary-layer solutions for various 
values of m (from m = -0.09 to m = 1) have been tabulated in detail. 
It has been assumed in these solutions that the viscosity, heat conduc- 
tivity, and specific-heat coefficients are proportional to certain powers 
of the absolute temperature . Thus, it is assumed that 


- ( T / T v) e ■ 

°p ./% - (*/*»)" 

The numerical values of o>, e, and k have been taken to be 


(27a) 


03 = 0.7 
e = 0.85 • 
k. = 0.19 


(27b) 


According to equations (27a) and (27b), the Prandtl number Pr£= \xc^jkj 
will vary very slightly according to the relation 


Pr/Pr w = (T/SV) 0 **" 6 = (T/Ty) 0 * 04 
It has been assumed in references 16 and 17 that 


(27c) 
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Pr v = 0.7 


(271) 


Since the pressure is assumed to be constant across the boundary- layer 
thickness, the ideal gas lav implies 

P/P-L = TjJT (28) 


Distribution of Normal Injection Velocity 

In order to obtain similar solutions of the boundary-layer partial 
differential equations for the cases represented by equation (25), a 
particular type of distribution of the normal injection mass flow along 
the vail must be assumed. In particular, it is assumed that f w is con- 
stant, where ^\r> which is dimensionless, is defined by 


f v 


2 

m + 1 


v v 



(29) 


From equation (29), in conjunction with equations (25), (27a), and (28), 
it follows that the distribution of normal injection mass flow along the 
wall will be 


<P(&) = 


p v v v 
Pb^ 





(50) 


(f w is negative for injection and positive for suction; l is a charac- 
teristic length in the direction of flow, not related to the boundary- 
layer length L in the stability equations). Thus the dimensionless 
injection mass-flow parameter 9 will be proportional to the constant -f w 
and will vary inversely as the (l - m)/2 power of the distance along the 
wall. 


For stagnation flows (m = 1), equation (30) shows that the normal 
injection mass flow at the wall implicit in the similar solutions treated 
here will be uniform. In that case, in fact. 


i 32 . 

Vl - - f w(®l/ T v) 


(31a) 


For flow over a flat plate (m = 0), on the other hand, 9 will diminish 
along the wall from an infinitely large value at the leading edge as shown 
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1—0} n /p 

Vo ' -(V^OVV)" 2 ' S' (»») 

Typical distributions of cp(|) are shown in figure 8. 

It may be noted that for a given (negative) value of f w and a 
given wall temperature ratio T w / Tp (< 1), the average value of 9 over 

the length l will be the same for all values of m. This follows from 
equation (30), according to which 9 will be independent of m (for 
m > -1). In particular, for m > -1, 


l-o> 



Equation (32) indicates the physical significance of the constant f w . 
Thus, -f-w is approximately equal to the average value along the wall 
of the injection mass -flow parameter 9. 

Equation (30), with f w constant, for the distribution of the nor- 
mal mass-flow parameter 9 has been obtained without consideration of 
any heat-balance equation between the fluid flowing over the wall and 
the cooling fluid injected through the porous wall. It can be shown, 
however, (cf . the section immediately following) that if Buch a heat 
balance is taken into account then the condition that the coolant tem- 
perature as well as the wall temperature be uniform yields the same type 
of distribution for 9(|) as does equation (30). 


Coolant Temperature 

In references 16 and 17, a heat-balance equation at the transpiration- 
cooled wall was not considered. Consequently, no relation is given there 
between the maintained uniform wall temperature and the mass-flow con- 
stant ^w on the one hand and the coolant temperature on the other hand. 
Such a relation, however, can be obtained without difficulty by first 
expressing the condition that the heat transferred by the hot gas to the 
wall is absorbed by the coolant. Thus, If air is the coolant, then 

(k aT/Sy) w = - li) (33) 

where T^ is the coolant temperature, which is assumed to be uniform along 
the wall. Solving equation (33) for T^ and making vise of equations (26) 
and (29) give the following egression for the coolant temperature ratio: 
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% 

*1 






where « w * = (dfl/dri) w and d = (T - T w )/(Tp - T w ). 


(34) 


Prom, equation (34) , the (implicit) value of Tp/Tp corresponding 
to a pair of values of T W/ /Tp(< l) and f w (< 0) can be determined for 
any case. 7 Corresponding values of the coolant temperature ratio are 
included in figure 8. 


Application of Tabulated Solutions to Stability Calculations 

In references 16 and 17, for a given value of m, the solutions of 
the ordinary differential equations for the velocity and temperature pro- 
files are tabulated explicitly for values of q throughout the boundary- 
layer thickness. Moreover, first, second, and third derivatives of these 
solutions are also tabulated, and from them the derivatives of the veloc- 
ity and temperature profiles required in the stability calculations can 
be easily obtained. For this reason the solutions tabulated in refer- 
ences 16 and 17 are especially well suited for the stability calculations 
to be performed in the present investigation. 

In applying equations (19) and (20) to the boundary-layer profiles 
tabulated in references 1 6 and 17, the boundary-layer length L has been 
taken to be the momentum thickness 5p. Application of equations (19) 
or (20) then yields directly the minimum critical Reynolds numbers based 
on the momentum thickness and the local velocity up. Observing that, 

for the boundary layers analyzed here, the momentum thickness can be 
expressed in the form 


S 


i 



(35) 


where K is independent of x (but depends on m, and f w ) , 

the minimum critical Reynolds numbers based on the free -stream length l 
and the reference conditions (i.e., at point b) can be obtained. Details 
of the explicit equations used in the calculations are given in appendix B. 


Tfor several of the cases treated in references 16 and 17 , negative 
values of Tp/Tp are obtained. This simply indicates that in order to 

maintain the given wall temperature ratio T w /Tp, a higher rate of mass- 
flow injection (i.e., a higher value of ”^w) would have to be used in 
actual practice. 
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STABILITY OF STAGNATION FLOWS AND FLOW 
OVER FLAT PIATE 


In order to investigate the stability of boundary- layer flows with 
normal fluid injection in a pressure gradient, the minimum critical 
Reynolds numbers have been calculated for two-dimensional flows in the 
immediate vicinity of a stagnation point (eq. (25) with m = 1) . This 
is a relatively simple physically significant flow with a favorable 
(negative) pressure gradient. In order to determine the effect of the 
pressure gradient, the minimum critical Reynolds numbers have also been 
calculated for flow in a zero pressure gradient or for flow over a flat 
plate (eq. (25) with m = 0) . For both cases, the tabulated solutions 
in references 16 and 17, as described in the section entitled "Mean-Flow 
or Steady-State Solutions," have been used. For each value of m, nine 
different flows have been tabulated in references 16 and 17, corresponding 
to three different wall temperature ratios (Tv/Ti = 1, l/2, and l/4) and 
three different values of the mass-flow injection parameter f v 
(f w = 0, -l/2, and -1) . The main results of the stability calculations 
are given In tables I and II and in figures 7 and 9 to 13 . 


Effect of Cooling of Wall and of Normal Mass-Flow Injection 

In figure 9, the minimum critical Reynolds numbers (^61) based 

min 

on the momentum thickness are shown as functions of the wall temperature 
ratio T W/ /T2_ for various fixed values of the normal injection mass-flow 
parameter f w .® The results for both stagnation flows and flow over a 
flat plate are shown in. this figure. For a given value of -f w it is 

seen that, as might be anticipated, the minimum critical Reynolds number 
increases, that is, the flow becomes more stable as the wall Is cooled 
or as T w yT]_ diminishes. This stabilizing effect of cooling of the wall 
(for a fixed value of -f w ) is seen here to be particularly pronounced 
for the flow over a flat plate at zero or moderate (^w = -l/2) rates of 

normal mass flow. It may also be observed from the nature of the curves 
in figure 9 that the stabilizing effect of cooling of the wall with a 
given value of f w is enhanced at the lower wall temperatures. 

®The square root of the minimum critical Reynolds number J 

Is plotted here on semilogarithmic paper merely for convenience in showing 
all of the results on one scale. Each curve for a' fixed value of f w is 
based on only three calculated points. A smooth curve through these pointe 
has nevertheless been drawn in order to show clearly the trends. See 
also table I. -• 
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In figure 10 , the Tniirfimnn critical Reynolds number for both stag- 
nation and flat -plate flows is shown as a function of the normal mass- 
flow parameter (cf. eq. ( 30 )) 

[2/(m + l)]<p| (l " m)/2 = -f w (^wAl)“ (l " C3)/2 

for various fixed values of the wall temperature ratio I ■ j T^'. Here the 

separate destabilizing effect of the normal mass flow can be clearly- 
seen. Ibis effect is seen to be especially pronounced for flow over a 
flat plate with a cooled wall (T v JT-]_ = l/2 and l/4). 

The net effect of simultaneous cooling of the wall and normal fluid 
injection can be clearly seen from figures 9 and 10 by comparing the 
results with those for the case of an impermeable uncooled wall 
(f v = 0, T w /Ti = l) . The latter case, for the stagnation and flat-plate 

flows, is shown by the corresponding horizontal lines in figures 9 and 10. 
For each type of flow (i.e., stagnation or flat-plate flow) a net sta- 
bilizing effect of transpiration cooling is indicated by those pairs of 
values of f w and T w /Tp or tp and T w /T]_ on the curves in figures 9 
of 10 which lie above the horizontal line corresponding to f w ~ 0 and 
T w /Ti = 1. Thus, particularly from figure 10, it is seen that for the 
stagnation flows the net effect of transpiration cooling (with T w /Tj_ < l) 

is stabilizing in all of the cases calculated herein, except for the case 
f w = -1 and T w /Tp = l/2, where the net effect is slightly destabilizing. 
Such is also the case for the flow over a flat plate (cf. fig. 10), except 
that, at the low temperature ratio T w /Tj_ = l/4, the high mass-flow injec- 
tion rate (f w = -l) practically cancels, or more than cancels, the sepa- 
rate stabilizing effect of the cooling of the wall. Flows over a highly 
cooled wall (e.g., = l/4) with low rates of normal mass flow (e.g., 

-f w ts 1 / 2 ) will, of course, be particularly stable and will have relatively 
high minimum critical Reynolds numbers. 

From figure 11 it is seen that the results for the minimum critical 
Reynolds numbers based on free-stream or reference conditions are essen- 
tially similar to the above results (cf . fig. 9) based on the momentum 
thickness and local conditions, provided that the free-stream Reynolds 

numbers are platted in the form (R^)^ 1°"" (see also eq. (Bll)). 


Effect of Pressure Gradient 

The effect of the favorable pressure gradient inherent in the stag- 
nation flows studied herein can be determined by comparing the minimum 
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critical Reynolds numbers for the stagnation flows with the corresponding 
minimum critical Reynolds numbers for flow over a flat plate. From fig- 
ure 9 or figure 10 it is thus seen that in most of the cases calculated 
herein the stagnation flow, for the same pair of values of ^w and T w /Tp, 
is more stable than the flow over a flat plate. The stabilizing effect of 
the favorable pressure gradient here is particularly great at the high 
injection rates (f v = -l) regardless of the wall temperature ratio V T i> 
or over a relatively uncooled wall (T w /T]_ « l) regardless of the injection 
rate, for the range considered herein. 

Most of the overall results thus far discussed might have been pre- 
dicted qualitatively on physical grounds. Thus, the separate stabilizing 
effect of cooling of the wall, the separate destabilizing effect of normal 
mass-floV injection, and the stabilizing effect, in most of the cases, of 
the favorable pressure gradient appear in accord with physical expecta- 
tions. A more detailed study of the results in table I and figures 9 
and 10, however, reveals several interesting features which might, per- 
haps, have been rather difficult to anticipate intuitively. 

From figure 9 , it is seen that except for the case of large rates of 
mass-flow injection (f w = -1) cooling of the wall with a fixed mass-flow 
parameter f w has a considerably greater stabilizing effect on the flow 
over a flat plate than on stagnation flows. In fact, figure 9 indicates 
that, at low wadi, temperature ratios (T w /Tp = l/4), (Rg.jJ for stag- 
nation flows will, as a result, not be higher than that for the corre- 
sponding flow over a flat plate. Thus, the stabilizing effect of the 
favorable pressure gradient here is greatly diminished and even elimi- 
nated at low wall temperature ratios and simultaneous zero or low 
(-f w 1a 1 / 2 ) rates of normal mass flow. These results may, to same extent, 

be explained by a comparison of the velocity profiles for the various 
cases. Velocity profiles are shown in figures 12(a) and 12(b), together 
with the corresponding values of c Q (cf. also tables Il(a) and 11(b)) 
on which the minimum critical Reynolds numbers strongly depend. Fig- 
ures 12(a) and 12(b) indicate that the velocity profiles for the stagna- 
tion flows are all fairly similar near the origin (and elsewhere) and 
that all have concave curvature, that is, negative values of cAi/dy2. 
Consequently, it may be expected that the minimum critical Reynolds num- 
bers will all be of the same order of magnitude and will vary relatively 
slowly with wall temperature and normal mass flow. For flow over a flat 
plate, on the other hand, figures 12(a) and 12(b) reveal that the shape 
of the velocity profiles changes appreciably as the wall is cooled from 
T w /Ti = 1 to T w /Tp = l/4 with a fixed value of f w . This can be seen 

especially clearly for the case of f v = -l/2. In this case, figure 12(a) 
shows that for an uncooled wall (T w /T]_ = 1) the velocity profile near the 
origin has an appreciable convex curvature , that is, positive values of 
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d'ni/dy^. The latter type of curvature tends to make values of J(c) 
algebraically small near the origin (cf. eq. ( 5 )) and hence tends to 
make values of c Q large. For f w = -l/2 with a highly cooled wall 
(T v /Tl = l/4 ) , however, it is seen from figure 12(b) that the velocity 
profile now has a concave curvature, that is, negative values of S^u/dy^ 
with consequently algebraically higher values of j(c) near the origin 
and a lower value of c 0 . The relatively large effect of lowered wall 
temperature on the stability of the flat-plate flows at low rates of nor- 
mal fluid injection can be further seen from figure 7* Here, for the 
case of zero normal mass flow, the curves J(c) and J m -; n (c) and their 
intersection ^determining 7 ] Cq and hence c Q ^ are shown for flat -plate 

and stagnation flows over an uncooled wall (T^/Tp = 1) and over a highly 
cooled wall (T w /Tp = l/4) . For the flat-plate flow, the considerable 
difference in the J(c) curves is seen to account for the much higher 
value of e Q ^and, hence, for the much lower value of ^ in the 

case T wAl = 1 than in the case T w /Tp = l/4. For the stagnation flow, 
on the other hand, the J(c) curves for T w /Tp = 1 and T w /Tp = l/4 
are seen to remain similar. 

It has already been pointed out that at zero or low rates of mass- 
flow injection, the stagnation flows, despite their favorable pressure 
gradient, will not be more stable and may even be less stable than flow 
over a flat plate when the ratio of wall to free-stream static temperature 
is very low. This perhaps unexpected result Is contrary to the trend pre- 
dicted by the results of all previous stability calculations involving 
a pressure gradient which have invariably indicated in one manner or 
another the stabilizing effect of a favorable pressure gradient and the 
destabilizing effect of an adverse pressure gradient. The present result, 
however, does not contradict any of these previous results since It g 
with a particular flow not included in the previous results. Thus, 
(assuming the correctness of the tabulated solutions In refs. 16 and 17 
and of the stability criteria used herein) the present result simply 
shows that there exists at least one type of flow with a favorable pres- 
sure gradient, namely a two-dimensional low-speed (Mp = 0 ) stagnation 
flow over a highly cooled impermeable or only slightly porous wall, which 
has a lower minimum, critical Reynolds number based on momentum thickness 
than the corresponding flow without a pressure gradient. A partial expla- 
nation of this result can be given by comparing velocity profiles and 
noting from figure 12(b), for the case of T w /Tp = l/4, that near the 

origin the velocity profiles for the flat -plate flow behave similarly to 
those for the stagnation flow when f w = 0 or -l /2 (in contrast, e.g., 

with the considerable change of the flat -plate profile when f w = - 1 , 

or with the different appearance of the flat -plate velocity profiles from 
that of the stagnation profiles when the wall is uncooled (cf . figs. 12 (a) 
and 7 ) • 
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From a practical viewpoint, the fact that the stability limit of the 
stagnation flow is apparently lower than that of the flat -plate flow at 
low or zero rates of injection over a highly cooled wall does not appear 

very serious . First of all, the difference in the values of ( Rg., ) 

' 1/ min 

is not very great. 9 Moreover, the minimum critical Reynolds numbers in 
both cases are so high that both types of flow may, for practical purposes, 
be regarded as completely stable. The present result nevertheless appears 
to be of considerable theoretical interest while it may also be of practi- 
cal interest if other cases of an apparent destabilizing nature of a favor- 
able pressure gradient exist. It would therefore appear worthwhile to 
investigate this problem further. For example, a recalculation of the 
present results based on mean-flow solutions obtained independently of 
those in references 16 and 17 may be worthwhile. Moreover, the extension 
of the present results for stagnation and flat-plate flows to still more 
highly cooled walls (T w /Tp < 1/4) would be of interest to determine whether 
the trends indicated in figure 9 continue as T w /Tp is decreased beyond 
l/4. Finally, an analysis of the stability of low-speed flows with pres- 
sure gradients other than stagnation flows, such as flows characterized 
by u x /u b = 1 ± a| (where a is a positive constant), over a highly 

cooled wall might reveal whether the result obtained herein is peculiar 
only to special types of flows such as stagnation flows. A calculation 
of the stability of three-dimensional stagnation flows (corresponding to 
m = l/3) would also be of interest. 

From figure 10 (cf. also table I), it Is seen that, except for the 
case of the uncooled wall, normal fluid injection with a fixed wall tem- 
perature ratio T w /Ti has an appreciably greater destabilizing effect 

on flow over a flat plate than on a stagnation flow. This effect is 
especially pronounced for the highly cooled wall (T w /Tp = l/4) at the 
high rate of fluid injection (f w = -l). Once again this can be largely 

explained by referring to the velocity profiles in figure 12(b) 

(T w /Tp = 1/4) where, for example, the considerable change in the flat- 
plate velocity profile as f w changes from -l/2 to -1 can be clearly 
seen. The stagnation profiles, on the other hand, tend to remain similar 
as f w varies from 0 to -1, regardless of the wall temperature (cf . 
figs. 12(a) and 12(b)). It is interesting, however, to note that, for 
_ _ 

In this connection, it may also be noted that the minimum critical 
Reynolds number based on main-stream reference conditions (^b) m; {_ n will 

vary inversely as | 2 for stagnation flows but inversely as £ for flat- 
plate flows (cf. eq. (Bll)). Consequently, for the same value of (RSi 

(Rb)min lor stagnation flows would exceed that for flat-plate flows at 
| < 1. In fact, according to table I, for the interesting case T w /T b = l/4 
and f w = 0, (Rbijjj^n for the stagnation flows will exceed that for flow 
over a flat plate at | < 13,510/24,900 or 0.343. 
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an uncooled or only moderately cooled wall, the stability of the flat- 
plate flows, in contrast with that of the stagnation flows, tends to 
become relatively less and less sensitive to changes in the normal mass- 
flow injection when the rate of the latter is high (ef. fig. 10). (This 
appears, in part, to be due mathematically to the large negative values 
of A for flat-plate flows with high injection rates, which tend to pre- 
vent (Rg. ) from diminishing sharply (eqs. (19) or (20)). 
v -‘-'min 


Effect of Coolant Temperature 

For a fixed rate of mass-flow injection cp(|), which corresponds 
approximately to a fixed value of f w (for a given value of m and a 

fixed value of <p(|), depends slightly on the wall temperature 

ratio T w /Ti (eqs. (30) to (31b)), lowering the coolant temperature Tp 
will lower the wall temperature T w . This, as has Just been seen, will 
stabilize the boundary layer by increasing the minimum critical Reynolds 
number. An example of this effect is shown in figure 13, where, for a 
fixed value of -1), the minimum critical Reynolds number is given 

as a function of the coolant temperature ratio Tp/Tq for both flat- 
plate and stagnation flows. 

It might also be of interest to calculate the stability of the 
boundary layer for a fixed coolant temperature ratio Tq/Tq and a varying 
rate of mass-flow injection. This has been carried out in reference 9 for 
(low-speed) flow over a flat plate. It has been shown there that for a 
fixed value of Tq/Tq there will be a value of the injection rate (essen- 
tially f w ) for which the minimum critical Reynolds number will be a maxi- 
mum. A similar calculation would be of interest for stagnation flows. 
However, boundary- layer solutions in addition to those tabulated in refer- 
ences 16 and 17 would be required for this purpose. One possibility might 
be to apply consistently the approximate solutions in reference 30 . 


Comparison of Stability Criteria 

The stability of all of the flows studied herein has been calculated 
by means of equations (8) and ( 19 ) • In view, however, of the quite recent 
appearance of reference 23 , a majority of the cases have been recalculated, 
for comparison, on the basis of equations (10) and (20), which are given 
in reference 23 as a modification, for greater accuracy, of equations (8) 
and ( 19 )* It has already been seen that these two sets of equations are, 
on the whole, essentially similar, especially for low-speed flows and 
moderate temperature changes within the boundary layer (cf. the section 
entitled "Stability Criteria"). The results of the calculations based on 
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equations (10) and (20) are shown in tables II (a) and II (b) and in fig- 
ure 10. Figure 10 indicates that, although some quantitative differences 
occur, the conclusions based on the results of equations (8) and ( 19 ) 
remain on the whole unchanged. It is, furthermore, significant to note, 
from tables Il(a) and Il(b), that the values of A calculated according 
to equation (10) are essentially similar to those based on equation (8) 
and that the calculated values of c Q consequently remain on the -whole 
quite similar. (When I'M « 1 , then the actual numerical value of A, 
whether it is positive or negative, will not have a great influence on 
the stability calculations . ) 

It is noteworthy that for the stagnation flowB studied here the 
absolute values of the profile stability parameter A are all much 
smaller than -unity (cf. table II (b)) regardless of the wall temperature 
ratio or of the mass-flow injection parameter f w . For the flow over a 
flat plate, however, relatively high negative values of A occur when 
the rate of mass-flow injection is high (cf. table 11(a)). This Is 
readily explained by a glance at the various velocity profiles in fig- 
ures 12(a) and 12(b) and a comparison of them with figure 6 (cf . also 
the section entitled "Stability Criteria"). The velocity profiles for 
the flat -plate flows at f v = -1 are thereby readily Been to be asso- 
ciated with high negative values of A. 


CONCLUSIONS 


The minimum critical Reynolds number as a function of the wall tem- 
perature and the normal mass -flow injection has been calculated for the 
laminar boundary layer over a transpiration-cooled wall maintained at a 
uniform temperature in the low-speed, but compressible, flow over a flat 
plate and in the immediate vicinity of a two-dimensional forward stagna- 
tion point. From these results, the following main conclusions have been 
drawn: 

(1) For a given rate of mass-flow injection, cooling of the wall 
stabilizes the flow by increasing the minimum critical Reynolds number . 

This effect Is enhanced at the lower wall temperatures. At a fixed wall 
temperature, an increase of normal mass-flow injection destabilizes the 
laminar boundary layer by decreasing the minimum critical Reynolds number. 
The net effect of simultaneous cooling of the wall and normal fluid injec- 
tion can be readily seen from the results obtained herein by comparing 
the minimum critical Reynolds numbers for this flow with those for flow 
over an uncooled impermeable wall. In most, though not all, of the partic- 
ular cases calculated herein, the net effect was found to be stabilizing. 

(2) In most of the cases, the favorable pressure gradient inherent 
in the stagnation flows has an appreciable stabilizing effect on the 
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laminar boundary layer. Ibis effect is especially pronounced at the 
high injection rates regardless of the vail temperature and over an 
uncooled or only slightly cooled wall regardless of the injection rate. 

(3) Except for the case of large rates of mass-flow injection, 
cooling of the wall, with a fixed rate of normal mass-flow injection, has 
a considerably greater stabilizing effect on the flow over a flat plate 
than on stagnation flows. At low ratios of wall temperature to free-stream 
temperature, in fact, the minimum critical Reynolds numbers for stagnation 
flows will consequently not be higher than that for the corresponding flow 
over a flat plate. Thus, the stabilizing effect of the favorable pressure 
gradient in stagnation flows is greatly diminished, and even eliminated 

or reversed, at low ratios of wall to free-stream temperature and simul- 
taneous zero or low rates of mass-flow injection. In view of this some- 
what unexpected result, it would be of interest to investigate the sta- 
bility of stagnation flows over still more highly cooled walls and to 
Investigate the stability of other types of flows with pressure gradients 
at low local speeds over a highly cooled wall. 

(4) Except for the case of an uncooled or very slightly cooled wall, 
normal fluid injection with a fixed wall temperature has a considerably 
greater destabilizing effect on flow over a flat plate than on a stagna- 
tion flow. This effect was found to be especially pronounced for the 
highly cooled wall at the high rates of fluid injection calculated herein. 

(5) For a fixed rate of mass-flow injection, lowering the coolant 
temperature will lower the wall temperature and hence will stabilize the 
boundary layer by increasing the minimum critical Reynolds number. 

(6) The velocity profiles for flow over a flat plate with high rates 
of mass-flow injection at the wall are of the type which lead to relatively 
large negative values of a stability parameter A. The approximate sta- 
bility criteria have been slightly modified herein for such cases. 


Polytechnic Institute 
Brooklyn, N. Y., 


of Brooklyn, 
November 1 6 , 


1955. 
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APPENDIX A 

COMPARISON OF EQUATIONS (8) AND (10) FOR BOUNDARY-LAYER 
PROFILE PARAMETER X(c) 


A qualitative comparison of equations (8) and (lO) for A(c) can 
be made by expanding the velocity profile w 3 u/u^ about the critical 

layer w = c (or y = y c ) . Thus, let 

v =* c + v c * (y - y c ) + - yc) + • • • (Al) 

where a prime Indicates 8/dy. To satisfy the condition u = 0 at 
y = 0, It is necessary, according to equation (Al), that 

w lf 

c - w c' y c + 2l" y c 2 ■ • • • ■ 0 (A2) 

Moreover, according to equation (Al), 



(u/ui) 



(A3) 


where £ = y c - y. Applying the binomial expansion in equation (A3) 
yields 


jT c ][ c - ( u / u i) dy = §(V) 1//2y c 3 ^ 2 - ^( w c , )’* 1//2w c ,,y c 5//2 + • • • 

(A4) 

From equation (Al), 

V = w c ’ - w c “y c + • • • (A5) 

Hence, when v/v 1 * 1 in the region 0 % y y c , equation (10) yields 
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A(c) 



Vc"y c + • • • )|2 ( .N 1 / 2 5 / 2 

■^572 i?(V) - 



-1/2 „ 5/2 


vJ'yo 





- 1 (A6) 


From equation (A2), 



1 + 


w c "y c 

2c 


(AT) 


Substitution for (w e 'y c /c) l n ^° equation (A6) yields, to the second 
power of y c . 





+ • • • 


(A8) 


With the use of equations (A5 ) and (A7), equation 
following expression for A(c), to the second power of 


(8) yields the 
y c : 


A(c) = -i I§_ y c 2 + . . . (A9) 

Equations (A8) and (A 9 ) show that, to second powers of y (and 
y c ) , equations (8) and (10) with v constant yield the same results for 

A(c) except for a slight difference in the numerical factors. Thus, it 
is seen that A and X depend on substantially similar properties of 
the velocity profiles and will both be of the same order of magnitude. 

If temperature changes within the boundary layer are taken into 
account and it is assumed that 


v/ v i = Vq/vj. + (v/vj *(y - y c ) + . . . (A10) 


then equation (A9) 


remains unchanged, while equation 


(A8) becomes 


A(c) 




+ w„ 




(All) 
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If v' > 0, as will be the case for low-speed flows with Vi < 1 > 

then the v terms in equation (10) will tend to make X algebraically 
less than X. This will tend to increase G^CX) (cf. fig. 5) or to 

decrease Jmin k ence to decrease c Q . 

In cases for which j X J « 1 (e.g., [Xj < 0.02),, it will be found 

that the stability criteria tend to be relatively insensitive to the 
actual numerical value of X. 
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APPENDIX B 

CALCULATION OF MINIMUM CRITICAL REYNOLDS NUMBER 


In the solutions of references 16 and 17* the functions f* f 1 * 
f" } f nt } d* fl 1 * and ■&" are tabulated* for any given case* as func- 
tions of tj } -where tj is defined by equation ( 26 ) in conjunction with 
equation (25). Primes denote here derivatives with respect to tj. The 
functions f and -9 are related to the velocity and temperature pro- 
files in accordance with the relations 


(Bl) 



From equations (Bl) and (B2) it follows that the derivatives of the 
velocity and temperature profiles can be calculated by means of the 
relations 


— (tj) = f 
U 1 


(n) 


' + *)• 



Since the expression for J(c)* as given by equation (5)* is inde- 
pendent of the boundary-layer length L* the differentiations indicated 
in equation (5) have been carried out with respect to tj. Similarly* 
in expressions (8) and (10) for A and A* which are also independent 
of L* y/L has been chosen to be tj. Thus* 


A(c) = (u/u 1 ) w ' % - 


(B4) 


1 
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Moreover, since in accordance with equations (27a) and (28) , 

. . l+o) 

V/Vi = (I/Tj.) 


(B5) 


equation (10) becomes 



l+O) 


( u / u l)w' 

c^ 7 ^" 



c - ( u / u l) 
/ , ■> 1 + 0 ) 
(T/Ti) 


drj - 1 


(b6 ) 


where 00 = 0.7* 

The expressions on either side of equations (19) and (20) do depend 
on the boundary-layer length L. For the present calculations, the 
length L has been chosen in these equations as the momentum thickness 
5 ± . From equations (26) and (35) it follows that 

t , = K(y/&i) (B7) 


Thus, equations (19) and (20), in conjunction with equation (B5), yield, 
respectively, the following expressions for the minimum critical Reynolds 
number based on the momentum thickness: 





i.^ 0 3( T / Tl )d«>[ 


2 Wti)(i - =o) 2 

\ v i ) 

min o 0 4 (u/u 1 ) Co ' (l + K Cq ) 5 

1 - 
— 

Vh - ^o ) 2 ] 1 ' 2 


(B8) 


\ E §±) 


min 


I u l 8 i \ 

\ V 1 / 


mm 


1 . 2 K 2 


[V T i) 

(2+0)), ^ 
( u / u 1 ) 

Ivl 1 ' 

==f 


(1 ^ *=o) 2 

1 - % 2 (i 

“ °o) 2 _ 

1/2 


(B9) 


In the. present calculations = 0. 


The minimum critical Reynolds numbers based on reference conditions 
(at point b) and the main-stream characteristic length l can be obtained 
from the results of equations (B8) and (B9) by using expression (35) for 
6 i . Thus, from equation (35) it follows that 



(BIO) 
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With | = x/z^ = v^ } and = u^l 131 , equation (BIO) in conjunction 
with equation (B5) yields 



From equations (b 8) or (B9) and (Bll) it follows that (Rh) . 

' 'min 

pendent of the momentum-thickness factor K. 


(Bll) 


is inde- 
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TABLE I.- MBUMtM CRITICAL REYNOIDS NUMBERS AS FUNCTIONS OF WALL 
TEMPERATURE AND RATE OP NORMAL MASS-FLOW INJECTION 


Stagnation flow. 

Flow over flat elate. 
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TABLE II— DETAILS OF SUABILITY CALCUIATIONS 
(a) Flow over a flat plate 



f 

Based on eqs. 

(8) and (19) 

Based on eqs. (10) and (20) 

Ti 

V 

mu 

MM 

( % i)min 

^0 

0 

O 


■j 

mm 

0.015 

0.415 

a 0.104 x lo3 

0.011 

0.413 

0.100 X k>3 


H9 

-.226 

.548 

.0225 

-.196 

•547 

.0232 

| 

mam 

-.703 

.620 

.0151 

-.663 

.622 

.0192 


mm 

0.035 

0.114 

12.5 




£ 


-.098 

.279 

.244 





mam 

-.702 

.459 

.0414 





mm 

msm 

40.0367 

378 

-0.009 

0.034 

486 


■9 

mmm 

4.0473 

85.0 

-.050 

.044 

108 


mam 

mm 

.319 

.100 

-.76 

.312 

.0665 


(b) Flow near a forward stagnation point 


Tw 

Ti 

f w 

Based on eqs. 

(8) and (19) 

Based on eqs. (10) and (20) 

n 

c o 

(" 1 U 

^0 

=0 

( ^®i)min 

n 

mm 

0.082 

0.210 

4.20 X lo5 

0.068 

0.207 

3-99 x 1£>5 


. 

.080 

.229 

2.43 

.055 

.223 

2.74 

mm 

mm 

.068 

.258 

1.36 

.058 

• 255 

1.56 


mm 

0.053 


22.4 




1 


.055 


8.62 





mm 

.046 

HgSSH 

2.49 





0 

0.049 

4o.o490 
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& Tbe accurate value of ( Rg., ) for this case, based on an exact calculation of neutral 
v 1/ min 

stability curve, without approximate criteria used herein, is 150 to 162 (refs. 20 and 18, respec- 
tively). Values based on particular simplified approximate criteria as given by eq. (21) range 
from 180 to 195 (refs. 12 and 20, respectively), which are tog high. Value obtained herein is 
evidently too low. These differences from exact value, however, are usually not regarded as very 
significant in stability calculations. It is important to note that in all of these calculations 
(i.e., refs. 12, IB, 20, and present report) for this case, value of c 0 is practically same. 

^Because of closeness of Cq to zero in these cases, determination of c 0 was based on 
expanding velocity and temperature profiles in a Taylor series about origin, using tabulated values 
of derivatives at origin and obtaining values of additional higher derivatives frets governing 
ordinary differential equations. ~ 

























































NACA TN 4037 


Figure 2.- Typical curve of G(z,A) versus z for a fixed value of A. 
Numbers indicate: 1, two values of z for chosen value of c sat- 
isfying equation (6); 2, one value of z for chosen value of c 
satisfying equation (6) ; and 3, no values of z for chosen value 
of c satisfying equation (6) . 
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Figure 4 .- Plot of z Q as function of A. 
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Figure J.- Determination of c Q for several different cases, t 
dots indicate location of critical layers. 
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Figure 8.- Distribution of normal mass-flow injection at wall 
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Figure 9*- Minimum critical Reynolds number, biased on momentum thickness, 
as function of wall temperature for various rates of normal fluid 
injection. 
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Figure 10.- Minimum critical Reynolds number Teased on momentum thickness 
as function of normal fluid injection rate for various wall temperatures. 
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Figure 11.- Minimum critical main -stream Reynolds number as function of 
wall temperature for various rates of normal fluid injection. 
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u/u, 


(b) T V /T L = 1/4. 


Figure 12.- Concluded. 







